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Abstract 

We study Fermi fields defined on tori in the presence of gauge 
backgrounds carrying non-trivial topology. We show that 2k dimen¬ 
sional field space can alternatively be described by fields over a k- 
dimensional space. This dual description is particularly natural when 

* Permanent Address. 
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the background is uniform and abelian. The reduction in number of 
dimensions carries over to the lattice. The lattice ultraviolet regu¬ 
larization induces an infrared regularization of the lower dimensional 
representations. We focus on A; = 1,2. 


1 Introduction 

Dirac fields defined over even dimensional tori can, for certain twists 
in the gauge background, be equivalently represented by fields de¬ 
fined over infinite flat spaces of half the number of dimensions. The 
two descriptions are related by a change of basis. The reduction in 
dimensionality comes from using Fourier space for half of the space- 
time variables. The integral Fourier momenta are used to string the 
remaining segments of continuous directions into infinite lines. This 
works for sections of non-trivial bundles. Although the two space- 
times are distinct in the topological sense, the fermionic index stays 
the same. For simple backgrounds, consisting of a constant abelian 
field strength, the field equations are local in either dimensionality. 

Extra interactions would be required to single out one of the alterna¬ 
tive descriptions as more natural. 

The setup we are looking at has a long history. It has been ana¬ 
lyzed before in various contexts [B El El m 0 El El- Early studies in 
lattice field theory were carried out in El. Probably the main objection 
to developing a lattice version was the absence of an exact fermionic 
“viewpoint” of background topology. Eermions react to topological 
features in the background in ways strongly correlated to their chiral 
properties and without a lattice version of chirality only an approxi¬ 
mate analysis was possible. This has changed with the domain wall 
formulation of lattice fermions [0] and, especially, with the overlap for¬ 
mulation |l()j . Among recent studies on the response of fermionic fields 
to topology in the overlap regularization we mention 

In this paper we focus mainly on the aspect of alternate description 
in terms of reduced dimensions. In the continuum, for our special 
backgrounds, locality is preserved and the equations of motion stay 
differential in either descriptionQ Our new contribution in this paper 
is to provide lattice analogues to the reduced dimension description. 

'^The description in terms of lower dimension fields is not restricted to Dirac fields: it 
would also be useful for analyzing the stability of the gauge backgrounds against small 
perturbations in the gauge fields mill 
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For example, we learn from this that the spectrum of lattice fermions 
moving on an Li x L 2 toroidal lattice in a uniform abelian gauge 
background depends only on the area L 1 L 2 . There is no distinction 
between the arbitrary Li x L 2 case and that on a degenerate lattice 
of dimensions L 1 L 2 x 1 . 

In sections 2 and 3 we introduce the continuum version of the 
representations in two and four dimensions respectively. The origin 
of the dimensional reduction is traced back to a natural Heisenberg 
algebra defined in this system. This algebra plays a central role in 
relating the developments of references na im to the the general 
area of Non-Commutative Field Theories, as reviewed, for example, 
in m- In section 4 we focus on the Dirac equation in constant abelian 
backgrounds, still in the continuum. Sections 5 and 6 deal with the 
lattice version in two and four dimensions respectively. We collect 
some numerical and asymptotic results in section 7. 

2 Two dimensions — continuum 

We consider a two dimensional torus of size li x I 2 . This can be 
viewed as the plane modulo the lattice generated by ei, 62 , chosen 
for simplicity as ei = (Zi,0) and 62 = ( 0 ,/ 2 )- The fields = 

'ip{xi,X 2 ) are complex fields satisfying certain periodicity properties. 
Our fields are coupled to a gauge field and are thus sections of a 1/(1) 
bundle on the 2-torus. The bundle can be twisted, the twist being 
characterized by the first Chern class ci. The first Chern number 
gives the topological charge: it measures the flux through the two 
torus in units of 27r (the charge e is set to 1 ). 

Periodicity of fermion fields has to be taken modulo gauge trans¬ 
formations: 

^|;{x + ei) = ^i{x)Tp{x) ( 1 ) 

'ip{x + e2) = D2(x)V’(x) (2) 

The U{1) fields are the transition functions made consistent by 

Di(x-I- 62 ) D 2 (x) = 112 ( 3 ; + ei) Di(x) (3) 

By gauge transformations the transition functions are brought to the 
form 

Di(x) = exp{^7^gx2//2} (4) 

D 2 (x) = exp{—z7r(7Xi//i} (5) 
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The consistency condition m forces q to take only integer values, q is 
the first Chern number (often referred to as winding number, in this 
set-up). 

Fields satisfying the boundary conditions (^I2I)make up a Hilbert 
space 7iq with scalar product 

{(pl'ip) = [ dx (j)*{x)'tp{x) (6) 

Jt'2 

The integration is over the torus and because of the integrand’s peri¬ 
odicity can be performed over any fundamental cell. 

TCq is left invariant by a ring of linear transformations generated 
by two unitary transformations Uj. 

Ui'ipix) = exp{2TTiXi/li}'ilj{x) (7) 

Ordinary translations map out of TCq because the transition func¬ 
tions depend on x. This can be compensated by an appropriate gauge 
transformation leading to modihed operations of translation by a*: 

Ti(ai)'I'(x) = exp{i'KqaiX2/A} ^{x + {ai,Q)) (8) 

7^(a2)^'(x) = exp{-z7rg'a2Xi/.A} T(x-I-(0, 02 )) (9) 

where A = I 1 I 2 is the area of the torus. The new translations along 
the two axes do not commute: 

71(01)7^(02) = exp{- 2 mqaia 2 /A} T2(a2)71(oi) (10) 

Elements of the non-commutative group of translations do not com¬ 
mute with the discrete abelian group of transformations generated by 
Uj. There exists only a finite subgroup of unitary transformations, 
generated by iF(^k(2) commutes with translations: 

ifW = ViT2{l2/q) (11) 

The boundary conditions imply = I and the operators KA) 

obey: 

^ exp{27rz/g} (13) 

Any one of the KA) generates a Zq group under which TCq decomposes 
into q orthogonal subspaces: 

7^9 “ ®n=l'^q,n (14) 


4 


( 15 ) 


Decomposing under K^‘^\ fermion fields in 'Hq^n satisfy; 

x) = exp{—27rz—} V'(x) 

K(i) 

maps these subspaces into each other, but the non-commutative 
translations leave the subspaces invariant. 

The generators of translations along each axis are given by: 

= di + mqeijXjjA (16) 

where €ij is the antisymmetric tensor with two indices(ei 2 = 1). The 
generators coincide with ordinary covariant derivatives in a constant 
magnetic field: 


Ai{x) = -nqeijXj/A (17) 

F = diA2 - 52^1 = ^ (18) 

The non-commutative translation group is made out of parallel trans¬ 
porters in this field, q is given by the total flux. The anti-self-adjoint 
generators of parallel transport do not commute; 

|I)r>,Dfl = -.F = ^ (19) 

Up to a rescaling, we have a Heisenberg algebra associated with one 
degree of freedom and we are led to a quantum mechanical view of 
the space as a Hilbert space of a particle in one-dimension. For this 
particle, the covariant derivatives play the role of position and momen¬ 
tum operators. For g > 0, we take P and Q ^ —iD^'^/F 

producing a traditional position space representation in terms of func¬ 
tions of a single real variable h{z). The space 7iq is identified with 
©gL 2 (R) with the D® action given by; 

hn{z) = iF z hn{z) (20) 

Dfhjz) = ( 21 ) 

The fields get mapped into the quantum mechanical position space 
basis as follows: Since the function 

X\X‘2 

ip{xi,X2) = expl-zTrg—V'(xi,X2) (22) 
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is periodic in the variable xi with period li we can expand it in a 
discrete Fourier series: 

ip{xi,X2) = ^exp{27r^^/c} hk{x2) (23) 

k 

The periodicity requirement in X 2 imposes a constraint on hk{x 2 )- 

hk+q{x 2 ) = hk{x 2 + h) (24) 

The constraint is solved by: 

hk{x 2 ) = hn{x 2 + kh/q) (25) 

with n = k mod q. The expression for the field becomes: 

1p{xi,X2) = explzTTQ'^^} X 

^ ^ exp{27rz|i(n + sg)} K{x 2 + (26) 

n=l sez 

The wave-function hn{z) transforms as required under and 
and the index n labels the subspace 'Hq,n to which the component 
belongs. 


3 Four dimensions — continuum 

We consider now fields defined over a four-dimensional torus of size 
I 1 XI 2 XI 3 XI 4 . If the gauge group is U (W) the bundles are classified by 
the instanton number Q and by the first Chern numbers associated 
to the individual planes. If the gauge group is SU(Nc) we only have 
the instanton number, but if all the fields are insensitive to the center 
Z{Nc) the true gauge group becomes SU{Nc)/Z{Nc) and there is a 
remnant of the first Chern numbers, now defined modulo Nc, labeling 
the so called twist sectors. Excluding the case W = 2 with zero twist 
and odd instanton number the topological character can be encoded 
in transition functions 12^(x) restricted to the 1/(1)'^'=“^ subgroup of 
su{Nc) mis]. The matter fields obey the boundary conditions: 

'tp{x= Qf_i{x)'tp{x) (27) 

The (not necessarily orthogonal) 4-vectors have lengths and 
span a non-degenerate four dimensional lattice A. The dual lattice 
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is generated by Denoting by the components of the 

invertible matrix R, ei^'^ = {R~^)ufj.- 

In the Abelian description we consider matter fields in the funda¬ 
mental representation with all Nc color components satisfying decou¬ 
pled boundary conditions: 

(x + (x) (x) (28) 

We again can pick transition matrices of the form: 

D““(x) = exp{z7rQ“^eJ,'^)xp} (29) 


The integers are antisymmetric in |U, i/ and can be assembled into 
the matrix Q°^. For U{Nc) the topological invariants are: 


a 

(30) 

Q = ^Pf(g“) 

a 

(31) 

Pf stands for Pfaffian. For SU{Nc), c^y = 0. 



We study the space of fields subject to the boundary conditions 
(UHl)-®- We deal with one of the components only and drop the color 
index; the number Q and the matrix Q°^ will become distinguishable 
only by context. We assume that Pf(Q) / 0. 

We again analyze classes of unitary operators that act in this space, 
starting with operations of multiplication by periodic functions. The 
generators of this discrete infinite abelian group are: 

= exp{27r^e[f^x^y} (32) 

A translation by 5, suitably modified to respect the boundary condi¬ 

tions, is given by: 

T{b)'ijj{x) = expi^FjPJb^Xy} V'(x -h b) (33) 

The antisymmetric matrix pjii) is given by: 

= 27rQ,.e(^)e(") (34) 

Translations only commute up to a phase: 

T(6)T(d) = expiiF^^Jbydf,} T{d)r{b) (35) 
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As in the two-dimensional case, we modify the ’s to achieve com¬ 
mutativity with translations. The generators of the resulting non- 
abelian group are given by {Q is invertible since by assumption its 
determinant does not vanish): 

k(p) ^ (36) 

They satisfy: 

= exp{-27r^(Q-^)^^} (37) 

The representations of this group, directly related to the Heisenberg 
group, have been studied in m (see m for a review). The irre¬ 
ducible representations have dimensions Pf((5). Our boundary con¬ 
ditions make the group not only finite, but also irreducible (elements 
commuting with all the generators are multiples of the identity) and 
the space of fields satisfying the boundary conditions (for a fixed color 
component) decomposes into a direct sum of Pf(Q) spaces. Trans¬ 
lations do not mix these subspaces. As we shall see later, there is 
exactly one zero-mode in each of these subspaces, giving a total of Q, 
as expected m- 

The infinitesimal generators of translations are given by 

= + (38) 

that is, covariant derivatives in a gauge field of uniform field strength: 

= (39) 

Before we derive a parametrisation of the space in terms of fields 
living in only two dimensions we exploit symmetries to simplify the 
problem. By transformations in SL{4:, Z) we bring to a canonical 
form with Q 12 = —Q 21 = Q and Q 03 = —Q 30 = q' and all other entries 
vanishing. The Pfaffian is equal to qq'. 

We now perform space-time coordinate transformations: 

^ = (40) 

In general, the transformation is non-orthogonal changing the metric 
tensor and the spectrum of the Dirac operator. However, the previous 
formulae remain valid, and the matrix F^h) is simpler now, with only 
Fm = -Ho =F' = and F,f = -F® =F = 2 vrq/(lil 2 ) 


as non-zero entries. As a result, the algebra of covariant derivatives 
decouples, becoming the Heisenberg algebra for a quantum particle in 
two dimensions and the problem factorises over the two planes 1 — 2 
and 0 — 3; so does the group generated by the operators . One 
can choose and as members of a complete set of commuting 

operators. Since = I the Hilbert space can be 

decomposed into subspaces labeled by integers 1 < n < 1 < n' < g'. 

Selecting —iD^f^/F and —iD^^/F' as position operators, we ex¬ 
press our fields as functions of two real variables, hn,n'{z, z'). The 
indices n,n' label the appropriate subspace. The covariant derivatives 
act as follows: 




z') 

= lFzhn{z,z') 

(41) 

Af /ini 

{z, 

z') 

= lF'zhn{z,z') 

(42) 


[z, 

z') 

II 

bi 

(43) 

Df /ini 

[z, 

z') 

= d^'hn{z,z) 

(44) 


The expression relating the initial four dimensional fields ip{x) and 
hn,n'{z,z') can be obtained from (f^. 

4 Uniform abelian fields 

We can use less dimensions for the fields, but, except for special back¬ 
grounds, the interactions will be non-local. If we are only concerned 
with fluctuations around one of the abelian backgrounds which played 
a role earlier, not only do the equations stay local, but, in addition, 
the Dirac operator becomes evidently exactly diagonalizable. As a re¬ 
sult, the famous relations between topology and fermionic zero modes 
is made very explicit. 

We begin in two dimensions with a gauge potential given by: 

Ai{x) = -Tiqeijixj - 2xf^)/A (45) 

The corresponding field is A = F 12 = 2T:q/{I 1 I 2 ). The coordinates 
x^f^ in Eq. specify the location where Polyakov loops around the 
torus take value 1 and represent non-trivial parameters of the gauge 
background. For x® = 0, the covariant derivative coincides with 
of section 2. In general, 

Di = di- iMx) = Ti(-xfV2(-xf (46) 
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Thus, the eigenvectors can always be obtained from those at x® = 0. 

The massless Dirac operator D = riD® + T 2 D^\ should have q 
zero modes. Let denote the Pauli matrices used to express the Dirac 
matrices 7 ^ in the chiral basis. The eigenvalue equation 


D 


( 0 ) 


D 


( 0 ) 




- 




^A 



(A) 

1 

(A) 

2 


(47) 


in the one-dimensional position space basis of section 2 is built out of: 

I)(0) ± W ^ ± Fz) = I (48) 

dz stands for the derivative with respect to 2 : and a is a standard 
harmonic oscillator annihilation operator: 


a 


a 


t 


1 

2 




1 

2 




(49) 


Dirac’s eigenvalue equation takes the form 



with representing the eigenmode wave function in the one dimen¬ 
sional basis. The index n takes q values and refers to the internal 
degree of freedom representing eigenstates of showing that spec¬ 
trum is g-fold degenerate. 

We express the expected positive chirality zero-mode in each sector 
in terms of \p), the eigenstate of the number operator a^a: 


| 0 ) 

0 


(51) 


The rest of the eigenstates are given by: 


where p = 1 , 2 .. 


hi^p) = ( 1 ^*) 

V2UIP-I) 


The corresponding eigenvalues are 


(52) 


( 53 ) 
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map into each other under chirality: 

In the original two-dimensional representation we obtain a well 
known form of the zero-mode. Starting from 



we get, with {x,y) G T^, 



(54) 


2g\V4 


= (ill) 


exp{^7^g 


X1X2 

~X' 


72=1 S^Z 


- 2 ^- htFxihin+sg) 


The Cn are the components of an arbitrary g-dimensional vector of 
unit modulus. 

The zero-mode is not spatially homogeneous in terms of either 
(xi,X 2 ) or z. In particular, for q = 1, the solution is proportional to a 
theta function. It vanishes at the point xi = /i/2,X2 = h/2: Eq. (l55|) 
gives 

4°) (xi,X2) = e*"/" ^ (-1)^ (56) 

The exponential is even under s —1 — s while (—1)^ is odd. The 
inhomogeneity of the zero-mode reflects the existence of the two extra 
physical parameters beyond the magnitude of the quan¬ 

tized magnetic flux, identifying where the Polyakov loops winding 
around the torus take value 1. At the coordinates of the zero of 
the two Polyakov loops take value -1. These coordinates could be used 
as parameters as well. 

According to Eqs. (15^ . (ISSl) the spectrum for g > 1 is the same 
as the spectrum of q independent fermions, each defined over a region 
of area A/q. This can be understood as follows; Dirac’s equation for 
q > I could be solved by using the q = I solution over the region 
0 < xi < h/q, 0 < X 2 < h (which carries one unit of magnetic 
flux) and extending it via the boundary condition '([’{xi + li/q,X 2 ) = 
exp[z7rx2//2]V’(2^i) 3 ^ 2 ) to the rest of the xi domain/i/g < xi < li. This 
gives the solution corresponding to the sector n = 0 , an eigenspace of 
Applying kA) one obtains the remaining q — 1 other solutions. 
All of them will satisfy the Dirac equation and the boundary condition 

V^(xi + h/q, X 2 ) = 3 . 3 ) ( 57 ) 
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Thus, the spectrum is that of q independent fermions defined over an 
area with a single unit of magnetic flux. 

Let us turn to the massive Dirac operator D + m. For g = 1 we 

get: 

( m \ ( \ 

Since in the presence of mass D + m is neither hermitian nor 
antihermitian it is better to consider the Hermitian Dirac operator 
H = t^D: 



m 

—i\/^a 



= E 


'hfy 

hf\ 


(59) 


H has one isolated eigenvalue at E = m with m-independent eigen- 
^ and an infinite set of paired eigenvalues, each pair 
labeled by p = 1, 2 ,3,.... 


function 


E^^P'> = ±^Jw? + 2pF, 


zi+p)z<^-p) = -1 (61) 


with 




_m jhh , / m^ hh , 7 

2 y pTT y 4 pvr 


(62) 


As in all other cases, the states and eigenvalues depend only on 
the area A = hh but not on the ratio h/h- One can view the one 
dimensional version as the limit when the ratio goes to zero or infinity 
while the area is kept fixed. 

When we vary m from positive values to negative ones isolated 
state crosses zero energy at m = 0. The paired eigenvalues never get 
to zero and are symmetric under m —m. The paired eigenstates 
depend on m, but the signs of z^^p') are m-independent. On the lattice 
there is full dependence on m, but the number of zero level crossings 
as m is varied stays the same as in continuum. 

Consider now four dimensions with gauge group SU{2) and fields 
in the fundamental, defined over an IqxIix12X torus. The boundary 
conditions are taken in the canonical form from the beginning; 


f ^ 

0 

0 

q'\ 

0 

0 

q 

0 

0 

-q 

0 

0 


0 

0 

0/ 


( 63 ) 
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± denotes the two color components. Note that the first Chern class 
invariants vanish. The topological charge is given hy Q = 2qq' with 
the factor 2 arises coming from the sum over color components. Odd 
topological charges cannot be obtained this way - as mentioned before, 
SU (2) is somewhat exceptional. The following constant field strength 
is compatible with the boundary conditions; 



f ° 

0 

0 

F'\ 

0 

0 

F 

0 

0 

-F 

0 

0 

\-T' 

0 

0 

0 / 


(64) 


where F = 2Trq/{lil2) and F' = 2TTq '/The upper color compo¬ 
nent coincides with the matrix F^^J of the previous section in canonical 
form. 

The fermion field is in the fundamental representation of SU{2) 
denoted by; 


T = 


• 0 + 

0 - 


(65) 


The index theorem tells us that there are Q = 2q(^ zero-modes of 
positive chirality more than zero modes of negative chirality. The 
eigenvalue equations decouple into color components. Because of the 
pseudoreality of SU{2) the upper color component determines the 
lower one by the symmetry 


T —> (66) 

act in spin and color space respectively while T* is the complex 
conjugate vector. The spectrum is doubly degenerate. We choose a 
chiral representation for the Dirac matrices; 


7m 


0 

(Tf, 0 J 


(67) 


75 = 


I 

0 



( 68 ) 


The matrices = (I, —it) are the Weyl matrices and are the 
adjoints. They satisfy; 




( 69 ) 
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is the ‘t Hooft symbol; = 6^^ and the are a basis of 
antiself-dual tensors. The eigenvalue problem is: 

Dipix) = = '^{x) = zAV’(x) (70) 


The symbol D = {D = ) denotes the positive (negative) 

chirality Weyl operator restricted to the upper color component. The 
covariant derivative coincides with the generator of translations 
studied in the last section, \ Choosing a new origin of coordinates 
gives the vector potential another form. The different choices are 
gauge inequivalent, as evidenced by the Polyakov loops, but are related 
by translations T(6). 

In position space representation D becomes {D = —D^): 

— _ / Do+iD^ i{Di-iD 2 )\ f iy/2F'a' - V^ a^ \ 

~ \ t{Di + zDa) Do-iDs ) ' V -V^a ) 

(71) 

a,a' are two independent annihilation operators, a is the same as in 
Eq. (jin|) . To get a' replace the 1 — 2 plane by the 0 — 3 plane and z 
by z'. 

A positive chirality eigenstate of the Dirac operator with eigenvalue 
zA, also is an eigenstate of the operator DD with eigenvalue — A^: 


DD = -2 


F'a'ta' + Fata 
0 


F'a'a'"^ + Faat ) 


(72) 


Let |p,y), with non-negative integers p and p', be simultaneous eigen¬ 
vectors of the number operators. There is a positive chirality zero¬ 
mode given by: 


/| 0 , 0 )\ 

0 

0 

\ 0 / 


(73) 


There is an identical zero-mode in each of the subspaces labeled by n- 
n' making up zero-modes for each color component and confirming 
the index theorem. Actually, all operators act in the same way in each 
subspace, so the entire spectrum is replicated Q times. 

The set of eigenvalues of the Dirac operator consists of: 


zA(^’P’P') = ±zx/2T>T2^ (74) 
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The corresponding eigenvectors are: 


1 

72 


I |P,P') \ 

_0 

^\/ Fp+F^ \P^P' ~ 

7>+fv ~ J 


(75) 


For p,p' > 0 there exist additional eigenstates with the same eigen¬ 
value: 


1 

72 




Fp+^F'p' \P^P' 

\ ^\lFp+F^ \P ~ ^^P') / 


(76) 


The spectrum and eigenstates of the massive hermitian Dirac op¬ 
erator H = 'y^{D + m) can be easily obtained in terms of those of the 
massless Dirac operator. The unchanged zero-modes become eigen¬ 
states of H with eigenvalue m. The massive Dirac operator mixes 
other eigenstates of the massless Dirac operator, ■0;^, with ysV’A- The 
eigenvalues and eigenvectors of H result from diagonalizing the corre¬ 
sponding 2x2 matrix giving eigenvalues equal to -|- A^. The 

eigenvectors are of the form T e*'^/^ 750 A); where d is the 

phase of the complex number m + ^A. 


5 Two dimensions — lattice 

We now turn to the lattice and examine a constant abelian background 
in two dimensions. We wish to determine how an ultra-violet cutoff 
interacts with the “dimensional reduction” we saw before. 

The fermion fields V’(ni,n 2 ) are complex functions of sites on an 
Li X L 2 square lattice. If we visualize the computer memory storing 
these functions, there is only a finite list of complex numbers and 
no mention of boundary conditions. No continuum features are lost 
because the link gauge variables are parallel transporters to a finite 
distance and therefore in the continuum limit include both the vector 
potential and the transition functions. The link variables thus keep 
track also of what we normally view as boundary conditions. To keep 
our lattice description formally close to the continuum formulation we 
choose to view our lattice fields slightly differently, as defined over an 
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infinite lattice but restricted by: 

ipin + = exp{nr'^qij^} ip{n) (77) 

j 

E^) is a vector of length Lj in the i direction, qij = qeij and the indices 
z,j take the values 1,2. 

The space of fields satisfying the boundary conditions (f77|) can be 
described in a way equivalent to the continuum one-dimensional quan¬ 
tum mechanical formulation. First, we gauge-transform the function 
V’(n): 

K{n) = exp{—V'(n) (78) 

L 1 T 2 

The new function, k, is periodic in ni with period Li and can be 
parametrised as: 

Li 

K{ni,n2)= ^ exp{^27^plnl/Ll}c(n2,Pl) (79) 

pi=i 

The momentum pi is defined modulo Li. Imposing the boundary 
condition in n 2 , we get: 

c{n 2 +L 2 ,pi) = c{n 2 ,pi+q) (80) 

Defining qi = gcd(g, Li), we can write pi = sq + n mod Li where n 
is an integer defined modulo qi and s is an integer defined modulo 
Li/qp. 

c{n2,pi) = h{n2 + L2S,n) (81) 

We have arrived at the following expression for ijj^n): 

V’(ni,n 2 ) = exp{7rzgf^} x 

Y.s=i^^ J2n=i exp{^27r(sg -h n)ni/Li} h{n2 + L 2 S,n) (82) 

For q=l, qi = 1, and our representation is given in terms of a function 
h defined on a one-dimensional lattice of length N = T 1 L 2 . For general 
g, it is still true that the space of fermions fields is Wdimensional, but 
now one has qi functions living on a periodic lattice of size N/qi. 

The lattice fields can be viewed as the collection of values 

of the continuum fields ip{x) (appearing in Section 2) at the lattice 
points X = na, where a is the lattice spacing and = L^o. Thus, 
the transformations acting on the space of lattice fields are given by 
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the restriction of the continuum transformations to the lattice points. 
In particular, we have a group of modified translations by vectors 
whose components are integer multiples of the lattice spacing a. The 
generators of the group are %: 

Ti'il;{n) = expjzTT^ 'ip{n + i) 

j 


(83) 


N equals L1L2 and z is a vector of unit length in the i direction. 
We obtained a non-commutative discrete subgroup of the continuum 
translation group: 


T1T2 = exp{- 27 rz^} T2T1 


(84) 


The one-dimensional representation is obtained in a basis that diago¬ 
nalizes Ti: 

Tih{k,n) = exp{27rz(^ -|- ^)} h{k,n) (85) 

N L\ 

The translation operator in the 2 direction acts as follows; 


T 2 h{k, n) = h{k + 1, n) 


( 86 ) 


If Ni 2 = gcd{N, q) ^ 1 the group generated by % acts reducibly 
on our space of lattice fields decomposing it into N 12 subspaces, each 
invariant under translations. This constitutes the lattice remnant of 
the group generated by in the continuum. These operators con¬ 
tain translations by li/q and can be extended to the lattice only if Lj 
is divisible by q. Otherwise, one obtains only a subgroup of the group 
found in the continuum. 

We pick lattice gauge fields having a constant plaquette value: 


Ui{n) = ai expjzTT^ 


N 


(87) 


The ai are complex numbers of unit modulus. The plaquettes are 
given by; 


Ui2{n) = exp{-z27r^} = C, ^ 


( 88 ) 


where Q = exp{z27r/A^}. The values of the plaquettes do not depend 
on cij, but the latter are physical parameters because they influence 
the Polyakov loops. The a^ are related to the continuum parameters 
X® (taken in units of the lattice spacing); 


— 2mq > . I 

ai = e 


.( 0 ) 


/N 


(89) 
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In the continuum the spectrum of the Dirac operators was independent 
of xf\ by unitary equivalence generated by appropriate translations. 
On the lattice translations are discrete and can only be used to restrict 
the ranges of the 0 < x® < 1. 

The lattice Dirac operators contain the covariant shift operators 
Ti-. 

Ti = aili (90) 

For the rest of the section we restrict ourselves to the <? = 1 case. 
The Ti operators take the following matrix form when acting on /i(/c)’s: 


0 


T 2 = 02 




1 

0 


0 

0 


0 ... 0 \ 

1 ... 0 


0 1 

0 0 / 


Ti — ai 


( 1 
0 
0 


0 0 
C 0 
0 


0 \ 

0 

0 


Vo 0 ... 0 


(91) 


(92) 


It is convenient to choose a different basis of Dirac matrices for the 
lattice analysis: 


7i = ^- 2,72 = T 3,73 = -^7172 = n (93) 

The Wilson Dirac operator D\y is 

Dw = m + Y,{r-V^) (94) 

where the unitary matrices are: 

V, = (95) 

In what follows we will fix the Wilson parameter r to 1. From now on 
periodicity in the “site” index n is implied by the convention that it 
is always viewed modulo N and this extends to operations involving 
site indices. For simplicity, let us start with x® = 0. 

The Wilson Dirac lattice Hamiltonian is 


Hw = tiDw 


/ — sin 0 

V 2 + m — cos 0 


Tt 
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2 + m — cos 0 — T 
sin 0 


(96) 




Here, © is an x diagonal matrix: 


G = diag[l,9i,02,...,ON-i] (97) 

with 

'^kn = h+l,n ( 98 ) 

Hw is a local one dimensional Hamiltonian defined on a discretised 
circle. The calculation of eigenvalues and eigenvectors requires fewer 
operations than in the basis corresponding to a two dimensional grid. 
Writing the eigenvector for energy E as 

^" = (t) <“> 

where n ranges between 0 and N — 1, we get the eigenvalue equations 

EUn = - sin Onttn + (2 + m - COS 6n)bn - bn+l (100) 

Ebn +1 = (2 + m - cos 9n+l)an+l - On + sin ^n+l^n+l (101) 


In matrix form we have 


2 + m - cos 9n+i sin 9n+i - f «n+l \ 
0 1 y V bn+l ) 

f I 0 

\—E — sin 9n 2 + m — cos 9n J \ b. 


Inverting the matrix on the left hand side we get 


The eigenvalue condition is that M = ■ ■ ■ Mq have an 

eigenvalue equal to unity. Observe that detM = 1; thus if M has 
one eigenvalue unity the other also must be unity and the eigenvalue 
condition can also be written as 


TrM{E) = 2 (104) 

When this condition is fulfilled M{E) has a single eigenvector, not 
two, because it is non-diagonalizable. The energies E are the roots of 
a polynomial of degree 2N. Typically these roots are all simple and to 
each corresponds a single eigenstate of H]v - there are no degeneracies. 
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We have ended up with an exact formula for the characteristic 
polynomial of Hw- 

det{E - Hw) = {-l)^KN[TrM{E) - 2] (105) 

The proportionality factor K]\f is calculated by comparing the terms 
E 2 N TrM{E) and in det(£; - Hw)- 


7V-1 

K]\f = (2 + m — cos = (106) 

n=0 


2 + m + ^y{2 + m)^ — 1 

N 

+ 

2+m—y/ (2+'m)2 —1 

N 

1 

2 

2 

2 N -1 


This formula is correct even when the square roots are imaginary, with 
the understanding that the second term on the right hand side is the 
complex conjugate of the first. It is derived by setting 2+m = ^{z+ \) 
and noting that the finite Laurent series in z defining iL^r is symmetric 
under z ^ ^ and vanishes when z is an root of unity. 

We can now easily re-introduce arbitrary x® and All that 

27rx^°^ 

happens is that 0 gets shifted to 0-and the characteristic 

polynomial is now proportional to 

TrM{E) — 2cos(27rx®) 


The proportionality constant now depends on 

Numerically, the finite precision of the machine will limit the use¬ 
fulness of a direct implementation of the characteristic polynomial as 
a means of locating the eigenvalues of Hw when N gets large. A bet¬ 
ter way exists for x^^^ = 0 which has extra symmetry. Ti and T2 are 
then isospectral and there exists an intertwining map, F]\f. (Actually, 
Ti,t|, T 2,r2 all have identical spectra.) F/v is the discrete Fourier 
transform: 


Fjk — 


y/N 




e N 


(107) 


F is symmetric and unitary. 


T 2 F = FTi (108) 

FT 2 = tIf (109) 
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F can be combined with a rotation in spinor space to produce a unitary 
matrix that commutes with Hw. 


U = 


1 - i 1 - i73 
^/2 ^/2 


( 110 ) 


One checks now that 


D\yU = UD\y (111) 

UHw = HwU (112) 

The simpler case of commutativity with U'^ is more obvious and can 
be checked directly; 


^ g -*0 


(113) 


leading to 


F‘^HwF‘^ = TiHwn (114) 

where we notationally ignore factors of unit matrices in tensor prod¬ 
ucts. 

The eigenvalues of U are l,i, —1, —i. Their degeneracies are com¬ 
puted from the ranks of the appropriate projectors. Since = 1 the 
projectors can be written in terms of 1, C/, U*, U'^. For the multiplici¬ 
ties we only need the traces of 1, F, F*, F"^. The trace of F^ is trivial 
because 


/I 

0 

0 

0 . 

. 0 


0 

0 

0 

0 . 

. 0 

1 

0 

■ • 0 

0 

0 . 

. 1 

0 

0 

• • 0 

1 

■ - 0 

. 0 

■ - 0 

Vo 

1 

0 

0 . 

. 0 

0/ 


The trace of F is not trivial and is given by a famous sum of Gauss: 

, Af-l 

Sg(N) = -^Y1 (116) 

SciN) is determined by its values at = 1, 2, 3,4 and by (nontrivial) 
mod 4 periodicity: <S'g(1) = 1, *S'g(2) = 0, S'g(3) = i, *S'g(4) = 1 -|- h 
The reflection U'^ = 73^^ = tiF‘^ reverses the order of sites with n = 0 
as a fixed point for all N and flips a and b (in our 7-matrix basis). 
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The origin of the ^/-symmetry is in the continuum. The harmonic 
oscillator case is unique in that all the eigenfunctions are eigenstates 
of the Fourier transform as a result of a discrete symmetry that ex¬ 
changes the coordinate with its conjugate momentum. The eigenval¬ 
ues of the Fourier transform are 1, i, —1, —i depending on the oscillator 
level mod 4. To get an operator that commutes with the Dirac oper¬ 
ator in our case the Fourier transform needs to be combined with the 
action of a rotation in spinor space depending on the chirality matrix 
73 . On the lattice one can decompose Hw into four blocks each of 
sizes (roughly) of N/2 x N/2. The exact sizes depend on N mod 4. 
The 17-symmetry is simpler to understand when also using the stan¬ 
dard representation in terms of fields on a two dimensional torus of 
equal sides. Then we have a symmetry under rotations by which 
corresponds to a switch between coordinate and momentum in the 
one-dimensional framework. The two dimensional interpretation as a 
discrete rotation explains why an action on spinorial indices is also 
required. As only the area is relevant, the requirement that the torus 
have equal sides can be dropped. In practice we shall exploit only 
corresponding to a rotation by tt, and the issue does not even arise. 
In the one dimensional basis the C/^ reflection symmetry contains the 
square of the Fourier transform combined with 73 in a simple manner. 

The 17^ symmetry allows us to factorize TrM(E) — 2. The trace 
contains the sum of the expectation of M{E) in two states. These 
contributions are individually high degree polynomials in E and vary 
over many orders of magnitude. The contributions nearly cancel at 
the eigenvalues and are very sensitive to E. The symmetry makes it 
possible to work in only one sector, where a single state contributes. 
As a result one deals with a single expectation value, which can be 
renormalized iteratively, in effect dividing out the quantity that has 
to vanish by a positive function of E that tames the large variability. 
This provides a numerically stable method for locating the eigenvalues 
down to machine accuracy. 

One needs to separate the case A^=even from A^=odd. For N = 2L, 
L-integer, the reflection has two fixed points: n = 0 and n = L. Thus, 
one can impose on ipQ and V'L to be eigenvectors of ri with identical 

= and let C- = ^ 

We need to have then 

CIMl-i{E)Ml-2{E) ... Mo{E)C+ = 0 (117) 

This gives a polynomial equation for E of degree 2L; the roots are all 


eigenvalue. Start from V’o = -^ 
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the energies of the even states. Similarly, for the odd states we have 
the equation: 

CiML-i{E)ML-2{E)... Mo{E)C- = 0 ( 118 ) 

Extending to = 2L — 1 requires imposing V’L-i = coordi¬ 

nated with zizTi'ipo = 'ipQ. This implies 


{1 - = 0 ( 119 ) 

if riV’o = V’o, and 

{1 + tiMl-i)^Pl-i = 0 ( 120 ) 

if TiV’o = -'4’0- 

Explicit evaluation of produces: 


Ml-1 


l-js+E)^ 


2+m+c 

-is + E) 


s + E 
2 + m + c 


( 121 ) 


where s = sin 2 l-i and c = cos 2 L-i ■ The eigenvalues of riMi_i are 
A = ±1, independently of E. The corresponding two eigenvectors do 
depend on E: 

4xiE) = , r/^ I F , w/n I I \l f ( A+g-Hs ) (122) 

The eigenvalue conditions are that (Mi_ 2 Mx ,_3 ... Mo).^o- be parallel 
to ia{E) for a = ±1. 

|+i(E')r2(Mi_2Mi_3 ... Mo)^+ = 0 (123) 

EYE)t2{Ml-2Ml-3 -.. Mo)e_ = 0 (124) 


They look slightly different from the even case because now the vectors 
^±1 are not orthogonal to each other. 

The above is useful for numerical searches for energy eigenvalues 
of Hw. The vectors obtained by the sequential action of are 
normalized at each step and this eliminates accuracy problems. N is 
limited now only because for very large N the spacings between some 
eigenvalues may be below machine accuracy. 

Looking at some examples the following patterns emerge: The odd 
and even states have energy sequences that separate (interlace with) 
each other. For even N the topological charge — rHw /is 
carried by the odd states. For odd N the topological charge is carried 
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by the odd and even states equally - each sector having an imbalance 
equal to unity. When m is taken from positive values to ~ —1 it 
is always an even state that does the zero crossing. There are no 
degeneracies. 

Since there are no degeneracies one can calculate directly the eigen¬ 
values A of the symmetry operator U for each one of the eigenstates. A 
simple way to do this is to evaluate F'tp at site 0 only. We distinguish 
the even and odd cases. 

If the state is even, 73 then 




iV —1 

yN 


and A can be ±1. If the state is odd, = —ip, then 




Af-l 


VNt 




n=0 


iXipo 


(125) 


(126) 


and i\ can be ± 1 . 

We now turn to the naive lattice Dirac operator. In the one¬ 
dimensional version we now have: 


Hn = 


-sine i(rt-r)A 
i(T-Tt) sine ) 


(127) 


If the Li are even there are two additional symmetries involving stag¬ 
gering operators along each direction: 


Siii{n) = (-l)”’V'(n) . (128) 

The Si commute among themselves, satisfy = 1 and also: 

SiT, = -T^Si . (129) 

The two anticommuting matrices Ilj = Si'ji, i = 1,2 commute with 
Hj\f. We diagonalize one of the Ilj’s, 112 for example. The eigenspaces 
of the two ±1 eigenvalues correspond to two species of staggered 
fermions. Going over to the one-dimensional description, the operator 
S 2 becomes the corresponding one dimensional staggering operator S: 

Snm = {-ir^nm (130) 

5^ = 1, 5 is hermitian and STS = —T. 
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We obtain a reduced Hamiltonian for naive fermions by introducing 
the following hermitian-unitary matrix: 

R = ^^I + ^^S2 (131) 

It brings H]\} into a simpler form: 

RHnR = -tsHnpi (132) 

where H]\fR = — 1 S 2 D 12 and 

Di2 = \{Ti- tI + ^^2 - iTI) (133) 

In the one dimensional representation we get: 

HNR = S{^{T-T^) + s\nQ) (134) 

The operator S 1 S 2 anticonmutes with Hjsfji, implying that the spec¬ 
trum of is real, symmetric around zero and doubly degenerate. 
We can cast the eigenvalue equation in the following form: 


Q-n+l \ 

bn+l ) 



(135) 


where 


(136) 


For E to be an eigenvalue, M{E) = M]\f{E) ■ ■ ■ Mi{E) must have an 
eigenvalue equal to 1. However, since the determinant of M{E) is 1, 
the condition becomes again: 


Tr{M{E)) = 2 


(137) 


An analysis similar to the one performed for the Wilson-Dirac operator 
would apply also here. 

Let us conclude this section by analyzing the existence and mul¬ 
tiplicity of zero modes. We restrict our attention to positive chirality 
modes. There always exists an equal number of negative chirality 
modes; this is a reflection of the famous phenomenon of fermion dou¬ 
bling. The equation for the zero mode is 

h^nh - hT-i = sin On (138) 
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This is a second order recursion relation determining all the elements 
of the series in terms of /ig^^ and . Our solutions must satisfy the 
periodicity condition = hn'^. It is necessary and sufficient that 

this condition be satisfied for two consecutive values of re. Since the 
problem is homogeneous this could lead to two conditions on the ratio 

/h^\ so existence is not guaranteed. To prove existence we make 
use of the fact that the solution, if it exists, satisfies hn'^ = For 
re = 1 this relation follows from the vanishing of 6q. For other values 
of re it can be proven by applying the iteration Eq. JISHI)- 

Most of our previous formulas have been derived for even N, so we 
will consider this case first. Then, ~ ^-Ar/2’ which is equivalent 
to periodicity at re = —N/ 2 . From the vanishing of sin 0^/2 we con¬ 
clude ~ ^7V/2-i-i’ which together with reflection around re = 0 

yields /i^]y-/2+i ~ ^7V/2+i (periodicity at — A ^/2 + 1 ). We conclude that 
the iteration procedure is automatically periodic for arbitrary values 
of and /i® and therefore there are two positive chirality and two 
negative chirality zero-modes. Their shape can be automatically de¬ 
termined by applying the iteration to any pair of initial values of /i® 
and Once one zero-mode has been determined the other can be 
found by applying the symmetry operator: 

( 139 ) 

In the odd N case, although staggering does not apply, the equa¬ 
tion for the zero-modes (USHl) and the reflection property (hn^ = h^l) 
still holds. This, however, is not enough to imply periodicity, but al¬ 
lows to reduce the two equations to one. This fixes the ratio 
A single zero-mode exists for each chirality. To conclude our analysis 
of zero-modes we mention that for the 7^ 0 case, the reflection 
property and the existence of zero-modes is lost. 

In conclusion, the one dimensional character of the lattice problem 
holds independently of which lattice fermions we use. Also, from the 
numerical point of view there are no hidden features and the problem 
truly has the complexity of a one dimensional system. 

6 Four dimensions — lattice 

In this section we look at an SU{ 2 ) gauge theory on a 4 dimensional 
lattice C of size Lg x Li x L2 x L3. We again take abelian transition 
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functions: 


ip{n + = expjzTT V tpin) (140) 

U ^ 

is an antisymmetric tensor whose only non-zero components are 
^703 = q' and qi 2 = q. The vector is a vector of length along 
the /U-direction. Following the two dimensional example we introdnce 
shift-operators which act as follows: 

T^V'+(n) = exp{^7r ^ q^lu^^} i’^{n + fi) (141) 

ly 

Changing the sign of the exponent on the left-hand side, one obtains 
the translation operator acting on the lower color component. The 
translation operators are just the covariant shift-operator in a back¬ 
ground lattice gauge field with the links given by: 

U^{n) = exp{z 7 r ^ (142) 

V 

These links produce a homogeneous plaquette field. 

The covariant shift operators define a finite gronp: 

= exp{-27r^-^} %% (143) 

For q = (^ = 1 the entire LoTiT 2 T 3 -dimensional space is irreducible 
under this group. A convenient basis is obtained by diagonalising a 
maximal set of commuting matrices; we choose Tj and Tq, leading to 
a two dimensional labeling of the basis given by the corresponding 
eigenvalnes of 7) and Tq. This is the lattice equivalent of the contin- 
unm two-dimensional representation of fonr dimensional lattice fields. 
Letting A^i 2 = T 1 L 2 and A'os = LqL^, the eigenvalnes of Ti are given 
by exp{z-^} and those of Tq are given by exp{z-^}, were m {m!) are 
integers defined modnlo A^i 2 (A^os)- An arbitrary vector in the space 
of lattice helds can be expanded in this basis with coefficients hm,m’ ■ 
From the commutation relations the action of the remaining operators 
T 2 and 73 can be derived by manipnlations from the previous section 
applied to each of the 1 — 2 and 0 — 3 planes. 

For the general q,q' case, the centralizer of the translation group 
might be non-trivial, leading to different subspaces invariant under 
translations. The snbspaces are labeled and mapped into each other 
by subgroups of the continuum group generated by the . If Li and 
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L 2 are divisible by q, and Lq and L 3 are divisible by we obtain qq' 
subspaces and labellings, just as in the continuum. On the opposite 
extreme, when the lengths are coprime with the associated integers 
q or q', the subgroup is trivial and we have no invariant subspace. In 
the general case we work with that subgroup of the continuum group 
generated by which maps lattice points into other lattice points. 

We set now q = q' = 1 and study the spectrum of four-dimensional 
lattice Dirac operators in the constant background field given by the 
links Eq. (CHI). We focus on the upper component, 'k+, because, as 
in the continuum, the lower component can be obtained by charge 
conjugation. 

The naive lattice Dirac operator is given by 

= ( 144 ) 


where we have defined the lattice Weyl operators in a way similar 
to the continuum. The non-zero spectrum of can be obtained 
by studying the eigenvalues and eigenvectors of the D^Dl acting on 
Weyl spinor lattice fields. The spectrum of this operator is the same 
as that of D^Dl. Unlike in the continuum, on the lattice the spectra 
are the same even including zero modes: for every positive chirality 
zero-mode there is a negative chirality one. 

In section 4 we related the spectrum of the 4-dimensional Dirac 
operator to that of the 2 -dimensional case and this relation extends to 
the lattice. 71^2 commute with and hence D^Dl becomes block- 
diagonal. 


DlDl 


DI2D12 + Djg 


Dq3 


0 


Di2D\2 + DosD, 


Ss, 


where 

Di 2 = l(T,-T^+^(T2-T^)) 

Dos = + 


(145) 


(146) 

(147) 


Di 2 and commute and our two dimensional analysis applies (see 
Eq. (11451) . Let A 12 and Aqs be the eigenvalues of DJ 2 D 12 and 
respectively. Then the spectrum of the naive lattice Dirac operator 
is given by Tz-v/AiJ^K^ and is two-fold degenerate. For the cor¬ 
responding naive Dirac Hamiltonian = 75 D 7 V, the spectrum is 



given by the same formula without the i factor. We emphasize that 
d\ 2 Di 2 = where is the staggered reduced Hamiltonian 

appearing in the last section above equation jni. Eigenvalues and 
eigenvectors can be readily constructed in terms of those of Hnr. 

The reduction of the 4-dimensional naive lattice Dirac operator 
to two dimensions does not extend to the Wilson-Dirac or the overlap 
operator. Here, all components are coupled to each other. A complete 
analysis, which will not be included here, can be done along the lines 
of the last section. 


7 Numerical and Asymptotic results 

In this section we will complement the results of the previous sections 
by providing some numerical and analytical results on the spectra of 
the naive lattice Dirac operator and the overlap Dirac operator |Tnll^ 
in a constant field strength gauge background for certain values of 
and topological charge. In particular, we will compare the eigenvectors 
of the lattice Dirac operators with their continuum counterparts. The 
numerical method used for the determination of the eigenvalues is an 
exact diagonalization procedure. This method, although precise and 
stable, is computationally costly, limiting what we can do. 

We begin by describing our two dimensional results. For the 
numerical determination we diagonalize the lattice Dirac operators 
(naive and overlap), with link variables given by: 

Ui{n) = exp{^7rg '^^^^’^^ } (148) 

^ 1-02 

U2{n) = expj-zTrg '^ (^i;^2) | 

L1L2 

where I{n 2 ,ni) = n 2 for 1 < ni < Li and /(n 2 , Li) = {Li + l)n 2 (and 
a similar relation for I' exchanging the 1-2 labels). The value of the 
plaquette is constant and equal to exp{—^27^■^}, with N = L 1 L 2 . The 
associated lattice Dirac operators (matrices) have dimension 2N. The 
cases of very small and very large N can be studied analytically. For 
N = 2 the naive Dirac operator vanishes. For N = A there are 
again 4 zero eigenvalues. In addition, there are two doubly-degenerate 
pairs at Ai.iy/2. For N = 6 there are three fourfold-degenerate eigen¬ 
values at 0 and ±z\/3/2. We have numerically diagonalized the naive 
lattice Dirac operator Djv for several values of q and Lj. In the ta¬ 
ble we display a selection of the eigenvalues for one typical case with 
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g = 1 and Li = L 2 = 24. We dropped the complex factor i from 
the eigenvalues and list only the positive branch on account of the 
spectral symmetry about zero. Since the continuum eigenvalues are 
given by y/2Fn, where F = 2'iTq/N is the magnetic held, we display 
the square of the eigenvalues divided by 2F^ a dimensionless quantity. 
We see that the degeneracy of the eigenvalues (indicated inside square 
brackets) is equal to 4. Only a factor of 2 is exact at all even values 
of N. The extra twofold degeneracy is approximate, but extremely 
accurate for large values of N. 

For large N the low lying eigenvalues are well described by: 

i\n = V^^{l + CFn + 0{F‘^)) (150) 

with C = —1/4. The dependence on N and q only enters through 
F and we are close to the continuum limit. C was determined by 
expanding D 12 to hrst subleading order in terms of continuum opera¬ 
tors: 

D\2 = F>i + iD2 + + 1 D 2 ) + ■ ■ ■ ( 151 ) 

0 

In the one-dimensional quantum mechanical representation we get: 

DI 2 D 12 = 2Fala “ af) +■■■ (152) 

The eigenvalues of —D\ 2 Di 2 are the square of the eigenvalues of Hat, 
and perturbation theory predicts the observed value for C in equa¬ 
tion (|150|) . 

We have also studied the eigenvalues and eigenvectors of the over¬ 
lap operator Dq- Dq is constructed in terms of the unitary operator 
n = T'iHwD{—M,r)/\HwD{—M^ir)\., where HwD{—M,r) is the Wil¬ 
son Dirac Hamiltonian with mass —M and Wilson parameter r. The 
eigenvalues of D are of the form ^ depen¬ 

dence explicitly indicated. Chirality implies symmetric spectrum and 
only the positive values will be shown. The continuum eigenvalues are 
approached by M6o{M,r) and results are presented for the square of 
this quantity divided by 2F. In terms of 6o{M,r) the eigenvalues of 
the overlap operator are: 

M(1 — cos 6o{M,r) + ism6o{M,r)) (153) 

This is the operator whose inverse gives the propagator for internal 
fermions lines. External lines should be described instead by the ex¬ 
actly chiral invariant operator Dj = 2iM^^ |in[l21| . The eigenvalues 
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of the latter are: 


i2M ian.{5o{M,r)/2) (154) 

Determining the eigenvalues for small values of N is somewhat 
more involved than for the naive operator, since now in addition we 
must keep track of the M and r dependence. For example for N = 2, 
q = I there is a region of values of these parameters for which D is equal 
to I and another one in which it is equal to —I. For 2r^—4Mr+M^ < 0 
(which includes the point r = M = 1) D has two eigenvalues equal to 1 
and two equal to —1. The table contains a selection of eigenvalues for 
(7 = 1, Li = L 2 = 24, r = 1 and M = 0.75,1,1.25. For large N, the 
behavior of the low lying spectrum can be understood analytically. 
The eigenvalues are again given by It™ where the constant C = 
C^overlap(-^;is. 

^overlap(Ikf, ^ “4 + (^^^) 

For r = 1 and M = 0.75,1,1.25 C = -11/108,1/12,37/300 respec¬ 
tively, which works well for the numbers in the table. 


Cont. 

Naive 

Ovp, M=1 

Ovp, M=0.75 

Ovp, M=1.25 

0 

0[4] 

0 

0 

0 

1 

0.9946 [4] 

1.0018 

0.9978 

1.0027 

2 

1.9783 [4] 

2.0071 

1.9910 

2.0108 

10 

9.4642 [4] 

10.167 

9.7723 

10.265 

20 

17.893 [4] 

20.615 

19.080 

21.040 

50 

37.371 [4] 

52.902 

44.316 

56.123 

100 

57.805 [4] 

106.45 

78.832 

121.57 

200 


199.11 

128.617 

259.15 

500 


350.07 

202.842 

523.46 


(fT^ is found by expanding the lattice operator in terms of the 
continuum one. We separate the Wilson-Dirac operator into additive 
terms as follows: 

Dwd = ~M + Dpf + rW (156) 

W stands for the Wilson term, which is of order a^, while the naive 
operator Hat is order a. a is the lattice spacing; equivalently, we 
set the lattice spacing to unity and use powers of instead. We 
write D = exp{—T/M}, where £ is an antihermitian matrix. In the 
continuum limit we get: 
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The above holds to order a^. In the continuum limit, W can be 
expressed in terms of the Dirac operator D: 

IT = + (158) 

Expressions (II .^71115SB are hold in any dimension and for arbitrary 
background. In two-dimensions the anticommutator of W and D is 
given by: 

WD + DW = -{D^ + t{D^F^,)Y) (159) 

The second term is proportional to the euclidean equations of motion 
and hence vanishes for any classical solution including the uniform 
fields we are considering. Hence, the behavior of the eigenvalues in the 
continuum limit can be extracted from our previously derived results 
for the naive lattice Dirac operator. In this way we obtain Eq. (Csni) 
with the C of Eq. dnsi)- The second and third term in Eq. (CSZI) 
are proportional to and therefore affect the eigenvalues without 
modifying the eigenvectors. 

For typical values of the parameters M, r (including r = M = 1) 
the lattice corrections are smaller for the overlap than for the naive 
Dirac operator. For example, for r = 1 and M = 2(1 — l/\/3) the 
leading 0{a?) correction to the eigenvalues vanishes. This also follows 
from the numerical data in the table, where the spectrum for r = 
M = 1 is found to be fairly close to continuum for the lowest several 
hundred eigenvalues. 

To conclude the description of the two-dimensional case we will 
show some results on the low-lying eigenvectors of the overlap (r = 
M = 1). In Figs 1 and 2 we display the shape of the positive chirality 
part of the ground state and the second excited state respectively. 
The data are obtained from the numerical diagonalization at Lj = 24. 
The resulting two-dimensional array of points is transformed to the 
one-dimensional representation and then compared to the continuum 
result, given by eigenstates of the Harmonic oscillator. 

We turn now to four dimensions. Our numerical work is done 
for SU{2) and for gauge fields whose only non-zero components are 
Fi 2 = —F 21 = Era and E 03 = —E 30 = EV 3 . The color components de¬ 
couple so we essentially are dealing with 17(1) systems. As explained 
in section 6 , the eigenvalue problem for the naive Dirac operator re¬ 
duces to the two-dimensional problem. Here we focus on the overlap 
operator, which does not simplify similarly. 
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Figure 1: The shape of the overlap (positive chirality) zero-mode in the 
one-dimensional representation compared to the continuum gaussian, for a 
2-dimensional constant held with q = 1 and Li = L 2 = 24. 
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Figure 2: The shape of the positive chirality component of the second excited 
state of the overlap compared with the continuum prediction, for g = 1 and 
Li = Z/2 = 24. 
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Our objective is to find the leading corrections to the low-lying 
spectrum when continuum is approached. Eqs. (HSU-dlSHl) are valid 
now, but Eq. (HSU) is replaced by: 

WD + DW = -{D^ + (160) 

where we have made use of the Bianchi identities. The second term 
vanishes for gauge configurations that solve the classical equations of 
motion. The last term was absent in two-dimensions. Eor a uniform 
held this term commutes with D. All the terms in £ anticommute 
with 75 and are antihermitian, so that in the chiral basis £ has the 
form; 

f=(" -^) ( 161 ) 

with X given by: 

X = Dn + - —)DDD 

^3M2 2M^ 

-^{F{Dq + ITsDs) + F’{lTiDi + ZT2T>2))t3 + ... (162) 

The eigenvalues of X^X are obtained by writing the covariant deriva¬ 
tives in terms of creation and annihilation operators and using per¬ 
turbation theory. With our choice of gamma matrices X^X is already 
diagonal in spin space. The 1 — 1 component of is readily ex¬ 

pressed in terms of the two number operators. The corresponding 
eigenvalue of X^X becomes: 

2Fni2 + 2F'no3 + F'n^sf 

- —(F^rei2 + T'^nos) (163) 

where ni 2 ,uo 3 are non-negative integers. The hrst two terms give the 
continuum result, the third one is of a type familiar from the two- 
dimensional case, while the fourth one is genuinely four-dimensional. 
The 2 — 2 spin-component of X^X gives another set of eigenvalues 
given by a formula identical to Eq. dMl) except that the last term 
has the opposite sign and ni 2 , no 3 are now integers strictly larger than 
zero. The sign change in the last term lifts the degeneracy of the 
continuum eigenvalues. The eigenvalues of £ are given by those of 
±^\/AtA. 

We have numerically determined the eigenvalues of the overlap op¬ 
erator in four-dimensions. Our exact diagonalization methods cannot 
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be carried out for very large sizes and a comparison with the previously 
derived formulas is difficult. 


8 Summary 

With one exception, any 17(1) or SU{N) gauge field on a two dimen¬ 
sional or four dimensional Euclidean torus can be smoothly deformed 
to an abelian background. When the background is topologically non¬ 
trivial the abelian field can be deformed to a non-vanishing uniform 
magnetic field. Any topological invariant is therefore captured by a 
representative with uniform abelian background. The Dirac equation 
in these backgrounds can be viewed as an equation over infinite space 
but in half the dimension. This naturally leads to an exact diagonal- 
ization in terms of harmonic oscillator wave functions. In the massless 
case one gets explicit formulas for the expected zero modes. Much of 
this translates to the lattice; for example the diagonalization of the 
Wilson Dirac operator and the associated overlap Dirac operator on 
a two dimensional torus simplifies to an equation on a one dimen¬ 
sional circle of length fixed by the area of the original torus. In these 
backgrounds the connection between fermions and topology becomes 
particularly transparent both in the continuum and on the lattice, 
once one employs overlap fermions there. 
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